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Abstract 

We present a phenomenological study of top squarks (^1,2) and bottom squarks 
(61^2) in the Minimal Supersymmetric Standard Model (MSSM) with complex pa- 
rameters At,Ab,iJ, and Mi. In particular we focus on the CP phase dependence of 
the branching ratios of ti^2 and 61^2 decays. We give the formulae of the two-body 
decay widths and present numerical results. We find that the effect of the phases 
on the ti^2 and 61^2 decays can be quite significant in a large region of the MSSM 
parameter space. This could have important implications for ti^2 and bi^2 searches 
and the MSSM parameter determination in future collider experiments. We have 
also estimated the accuracy expected in the determination of the parameters of ti 
and bi by a global fit of the measured masses, decay branching ratios and produc- 
tion cross sections at e~^e~ linear colliders with polarized beams. Analysing two 
scenarios, we find that the fundamental parameters apart from At and Af, can be 
determined with errors of 1 % to 2 %, assuming an integrated luminosity of 1 ab~^ 
and a sufficiently large c.m.s. energy to produce also the heavier ^2 and 62 states. 
The parameter At can be determined with an error of 2 - 3%, whereas the error on 
Ah is likely to be of the order of 50 %. 



1 Introduction 

Supersymmetry (SUSY) is one of the most attractive and best studied extensions of the 
standard model (SM) [Ij. With SUSY the hierarchy problem can be solved and the mass 
of the Higgs boson can be stabilized against radiative corrections. While this is certainly 
the main motivation, SUSY gives us the additional benefit of introducing potential new 
sources of CP violation [2 Ej- As the tiny amount of CP violation in the SM is not 
sufficient to explain the baryon asymmetry of the universe j3] , the systematic study of all 
implications of the complex SUSY parameters becomes absolutely necessary. 
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In the present paper we study the effects of complex SUSY parameters on the phe- 
nomenology of the scalar top quark and scalar bottom quark system. Analysing the 
properties of 3rd generation sfermions is particularly interesting, because of the effects 
of the large Yukawa couplings. Their lighter mass eigenstates may be among the light 
SUSY particles and they could be investigated at the Tevatron and at e^e~ linear colliders 
At LHC these states can be produced directly or in cascade decays of heavier 
SUSY particles [IHl^Iini- Analyses of the decays of the 3rd generation sfermions ^1^2, ^1,25 
fi 2 and Vr in the Minimal Supersymmetric Standard Model (MSSM) with real parame- 
ters have been performed in Refs. |I7j-jTn|. Phenomenological studies of production and 
decays of the 3rd generation sfermions at future e^e~ linear colliders, again in the real 
MSSM, have been made in Refs. 0^0- 

In the MSSM several SUSY breaking parameters and the higgsino mass parameter 
can be complex. In a complete phenomenological analysis of production and decays of 
third generation sfermions one has to take into account that the SUSY parameters Af, 
H and Mi [i = 1,2,3) are complex in general, where Aj is the trilinear scalar coupling 
parameter of the sfermion fi, and the Mi, M2 and M3 are the U(l), SU(2) and SU(3) 
gaugino mass parameters, respectively. This means that one has to study the effects of 
the phases of the parameters on all observables. 

An unambiguous signal for the CP phases would be provided by a measurement of a 
CP-odd observable. For example, in the case of sfermion decays a rate asymmetry [20] 
and triple product correlations [211 122] have been proposed as such observables. However, 
since it may be difficult to measure these CP-odd observables of the sfermions, CP-even 
observables like decay branching ratios may also be suitable to obtain informations about 
the SUSY CP phases. For example, the decay branching ratios of the Higgs bosons depend 
strongly on the complex phases of the t and b sectors ^23j-^25j, while those of the staus fi^2 
and r-sneutrino i>r can be quite sensitive to the phases of the stau and gaugino-higgsino 
sectors [20] • Also the Yukawa couplings of the third generation sfermions are sensitive 
to the SUSY phases at one- loop level [22] • Furthermore, explicit CP violation in the 
Higgs sector can be induced by i and b loops if the parameters At, Af, and /i are complex 
[231l2Hll2nilSni- It is found [23ll2nil2HllSI] that these CP phase effects could significantly 
influence the phenomenology of the Higgs boson sector. 

The experimental upper bounds on the electric dipole moments (EDM's) of electron, 
neutron and the ^^^Hg and ^os^l atoms may impose constraints on the size of the SUSY 
CP phases [H21ISn]- However, these constraints are highly model dependent. This means 
that the various SUSY CP phases need not necessarily be small. For instance, if we adopt 
the MSSM and assume a cancellation mechanism ^34], it turns out that the phase of fi is 
restricted as \(p^\ < vr/lO while the phases ipAf of the Af parameters are not constrained. 
On the other hand, the size of is not constrained by the EDM's in a model where 
the masses of the first and second generation sfermions are large (above the TeV scale) 
while the masses of the third generation sfermions are small (below 1 TeV) [201 • The 
restrictions on ip^ due to the electron EDM can also be circumvented if lepton flavour 
violating terms are present in the slepton sector [00] • Less restrictive constraints on the 
phases appear at two-loop level where 3rd generation sfermion loops can contribute to 
the EDM's [OZ]. 

In this article we focus on the influence of the CP violating SUSY phases on the 
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fermionic and bosonic two-body decay branching ratios of 3rd generation squarks ti^2 and 
bi,2- We use the MSSM as a general framework and we assume that the parameters At, 
Ab, /i and Mi are complex with phases ipAt, V^Aj,, f^i and respectively (taking M2,3 

real). We neglect flavor changing CP phases and assume that the squark mass matrices 
and trilinear scalar coupling parameters are flavor diagonal. We take into account the 
explicit CP violation in the Higgs sector. If the top squark and bottom squark decay 
branching ratios show an appreciable phase dependence, this would also affect the analyses 
of the various gluino cascade decays such as those in [1^]. In j2Hl we have published first 
results of our study. In the present paper we give the analytic expressions for the various 
decay widths for the complex parameters and study in detail the phase dependences of 
the branching ratios. We take into account the restrictions on the MSSM parameters 
from the experimental data on the rare decay 6 — > 57 [39 . Furthermore, we give a 
theoretical estimate of the precision expected for the determination of the complex top 
squark and bottom squark parameters by measuring suitable observables including the 
decay branching ratios in typical future collider experiments. 

In Section 2 we give the formulae necessary to calculate the U and bi two-body decay 
widths in the presence of CP phases. In Section 3 we present our numerical results. 
In Section 4 we give a theoretical estimate how precisely the complex top squark and 
bottom squark parameters can be determined at future collider experiments. We present 
our conclusions in Section 5. 



2 Squark masses, mixing and decay widths 
2.1 Masses and mixing in squark sector 

The left-right mixing of the top squarks and bottom squarks is described by a hermitian 
2x2 mass matrix which in the basis {qi, q^) reads 

M? Ml • ' - ' 



IRL IRR 



with 



= Ml + {II^-egSm^9w)cos2(3ml + ml, (2) 
^Lr = Ml, + e,sm^ewcos2(3ml + ml (3) 
ML = (M^J* = m,(A,-/i*(tan/5)-2^3.), (4) 

where rriq, Cq and I^i^ are the mass, electric charge and weak isospin of the quark q = 
b,t. Oy/ denotes the weak mixing angle, tan/? = V2/V1 with vi (^2) being the vacuum 
expectation value of the Higgs field (H^) and M^, = (M^) for g = 6 (t). Mg, 
M^, Mjj, Ah and At are the soft SUSY-breaking parameters of the top squark and bottom 
squark system. In the case of complex parameters /i and Ag the off-diagonal elements 
M? = (M| )* are also complex with the phase 



^q = arg [M|^J = arg Aq - /i*(tan/5)-2^3L . (5) 



3 



The mass eigenstates are 
with the g-mixing matrix 



, e*'^'? cos 6a sin Oa 
TV=\ ' . ' I , (7) 

— sin 6q e cos 6q 

cosg„- = ^= ' sin^, - 



,/|M| P + (m? -M| )2' " ./|M| P + (m? - M| )2 

and the mass eigenvalues 

I {ml + ML T ^iML-MlJ + mW) , ^91 < m,, . (9) 



"^?-1.2 



2.2 Fermionic decay widths of ti and 6j 

In the following we give the formulae necessary to calculate the two-body decay widths of 
ii and bi into charginos and neutralinos in the presence of the CP phases. The b — ii — x^ 
and t — bi — couplings are defined by 

>C99x+ =9i{t,^PR + ^.,PL)xth + 9h{e,^Pn + ^ P,) k + h-c. (10) 

with 

Pl = \{1-1,). Pr=\{1 + ^,), (11) 

and 

* ~ ^ ~ ^' ^ 

where g is the SU(2)^ gauge coupling and the 2x2 chargino mixing matrices U and V are 
defined in Eq. fEHjl . 

The q — (ji — Xk couplings {q = t,b) are defined by 

^99-x° = 9q{alkPR + ^1kPL)xlqi + ^-C- , (15) 

with 

2 2 

E(^-)*'^L> 4 = E(^-)*^^' (16) 



n=l n=l 



where 



K"'Rk/ \J Rk 



Lk 



Rk 
''Lk 



2^2 



^ ^^^^"^2 + 2 tan 9wNki ] , 



{h'RkY 



(1^ 



and 



fbk 



rb 
J Rk 

hlk 



{h'nkT 



(19) 



The 4x4 neutralino mixing matrix is defined in Eq. (I5()j) . 

The partial decay widths of qi {(ji = ti, hi) into fermionic final states then read 



r(g.^g' + Xfc) 



^2A2(m|,m2,,m2±) 

Hi 



i4r+i4r 



m^;, — m„, — m-± 

9« <? Xfe 



4Re(4*4)m,,m-± (20) 



and 



r(g.-^g + x°) 



^2A2(m| ,m^,m|o) 



IGvrm? 

,9 |2 , \rM |2 



X 



«r.r + ) ( < -ml- m|„ ) - 4Re(a?:6?,)m,m,o (21) 



with A(x, z) = + y"^ + — 2{xy + xz + yz). 
2.3 Bosonic decay widths of ti and hi 

Here we show the couplings relevant for the two-body decays of and hi into gauge and 



Higgs bosons. The qi — g' — W couplings are defined by 



c^^'w = -^9 ( d'^ h + ^r^. w, 



with 



bitj 



(22) 
(23) 



The qi — qj — Z interaction Lagrangian reads 



Cg^z = -^9Bf^Z^q]^^q, 



(24) 
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with 



/q r cos^ 9n — e„ sin^ 9i 



B^j — — 

COS tiw I 1 gjj^ 26, e'^^ sin^ - e, sin^ 



2 -"SL •^"^^>^g 

The gj — g' — i?^ couphngs are defined by 



with 



and 



^ ml tan /3 + m^^ cot /3 - sin 2/3 {\Ab\e'''^^b tan /3 + lAtje^^") A 



G 



\ 



mt {\At\e''^M cot + |//|e-^'^'') 



2mtmb 



sin 2/3 y 

For the couphngs of squarks to neutral Higgs bosons we have the Lagrangian 

^qqH ^-g C{qlHiqj ) qlHiqj {k,j ^1,2) 



with 



^ C{qiH,qL) C{qlH,qn) ^ 



C{qiH,qj) = 



\ C{qlH,qL) C{qlH,qR) ) 



where for q = t 



m 



mz fl 2 . 



mw sin /? cos ^vk \ 2 3 



sin 6w ) (cos l30u — sin l302i) , 



C{i^^H,iR) = 



7— 7;02i + ^ ^"^^ sin^ 9w{cosl30ii — sin /302i), 



m^y sin j3 3 cos 



(cos /3 1 At I e-^^^* +sinl3\^\e''P^) O, 



2mw sin /3 

|/x|e^^-Oi,-|At|e-^^-*02.)}, 



3i 



while for g = 6 



C{i^j,H,iL) = [C{i[HdR)]*, 



-Oii - "^f \\~\ si^^ 1 (cos/30ij - sin/302i) , 



■mwcos(5 cos^vK \2 3 



C(b^j,H,bR) = . Oi, - TT^^ sin^ OwicospOu - smP02i), (36) 

mw cos p 3 cos 0^ 



2mw cos p ^ 



3i 



-(l/ile^'^-Os.-IAIe-^^-^^^Oi,)}, (37) 



CiPj^Hih) = [C(b[H,bn)Y. (38) 

Here the 3x3 neutral Higgs mixing matrix O is defined in Eq. fl4(i|l . 

The partial decay widths for qi = ti, hi into bosonic final states are then of the following 
forms: 

m ^W^ + q ) = ^ , (39) 







167rm^m? 




\Bii\ 


2A5(m|,m|,m?J 






167rm|m?2 ' 






|2A5(m|,m|±,m?,) 



m ^H^ + %) = — ' 3 ' (41) 



3 Numerical results 

Before presenting numerical results, we briefly comment on the CP phase dependence of 
the (jiqj pair production cross sections. The reaction e~^e~ — * g^g^ (g^ = tj, bi) proceeds via 
7 and Z exchange in the s-channel. The Zqiqj couplings are defined in Eqs. ()24|) and ()25|) . 
The tree-level cross sections jHJ of the reactions e~^e~ qiq^ do not explicitly depend 
on the phases (p^ and ipA^- In the case of the reaction e+e^ — > g^gj, i = 1,2, the couplings 
Zg^gj are real. In e"*"e^ ^ gigj only Z exchange contributes and consequently the phase 
(fiq drops out in the matrix element squared. The tree-level cross sections depend only 
on the mass eigenvalues rriq^ ^ and on the mixing angle cos^ 9q. Therefore, they depend 
only implicitly on the phases via the cos(v?/i -l- v^aJ dependence of rriq^ ^ and 6q (Eqs. (jHI) 
and One-loop corrections to the cross sections have been calculated in for real 
parameters. In the energy range considered here they are of the order of 10%. For 
complex parameters they are expected to be of the same order of magnitude. Therefore 
we further expect that the direct influence of the phases on the cross sections as caused 
by one-loop corrections would be within a few percent. These phase effects on the cross 
sections would be much smaller than those on the tree-level decay widths studied in this 
paper. 
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In the following we will present numerical results for the phase dependences of the ti 
and bi partial decay widths and branching ratios. We calculate the partial decay widths in 
Born approximation according to the expressions given in the preceding section. In some 
cases the one-loop SUSY QCD corrections are important. The analyses of fHl iffl 112] 
suggest that a significant part of the one-loop SUSY QCD corrections to certain partial 
widths of ti and bi decays (where the bottom Yukawa coupling gYb is involved) can be 
incorporated by using an appropriately corrected bottom quark mass. In this spirit we 
calculate the tree-level widths of the ti and bi decays by using on-shell masses for the 
kinematic terms (such as a phase space factor) and by taking running t and b quark 
masses for the Yukawa couplings gYt^b- For definiteness we take m™'^(m^) = 150 GeV, 
^on-shell ^ GeV, ml'^^imz) = 3 GeV and m"^-^!^^" = 5 GeV. This approach leads to 
an "improved" Born approximation which takes into account an essential part of the one- 
loop SUSY QCD corrections to the ii and bi partial decay widths and predicts their phase 
dependences more accurately than the "naive" tree-level calculation. The inclusion of the 
full one-loop corrections to the partial decay widths of ti and bi is beyond the scope of the 
present paper. One-loop corrections to partial decay widths of ii and bi have been given 
in 1^ 133] for real MSSM parameters and are of the order of 10 %. We expect that for 
complex parameters they are of the same order of magnitude. In the calculation of the CP 
violating effects in the neutral Higgs sector we take the program FeynHiggs2.0.2 of [30j, 
which includes the full one-loop corrections to the mass eigenvalues and mixing matrix 
of the neutral Higgs bosons for complex parameters. For comparison we have also used 
the program cph.f of [2H]- We have found agreement between the results obtained with 
cph.f and the one-loop version of FeynHiggs2.0.2. There are small numerical differences 
between the results of cph.f and the two-loop version of FeynHiggs2.0.2. 

In the numerical analysis we impose the following conditions in order to fulfill the 
experimental and theoretical constraints: 

(i) m-± > 103 GeV, m^o > 50 GeV, m,- r > 100 GeV, m,- r > m^o, 

^ ' Xi ' Xi ' tijOi ' 11,01 Xi ' 

(ii) for incorporating the experimental bound on the mass of the lightest Higgs boson Hi 
we use Fig. 4 of [33], replacing mt by rriHi and sin^(/?— a) by (On cos/?+02i sin P)"^,^ 

(iii) 2.0 X 10~^ < B{b 57) < 4.5 x 10~^ |39| assuming the Kobayashi-Maskawa mixing 
also for the squark sector, 

(iv) Ap{i-b) < 0.0012 

(v) I < 3(M? + M? + ml), I < 3(M? + M? + mj) with 

Ic,/ C J-J 

m\ = (m^± -|- m| sin^ 9w) sin^ P — |m|, ml = (m^± + m| sin^ 9w) cos^ /? — 

Conditions (i) and (ii) are imposed to satisfy the experimental mass bounds from LEP 
133 133 EH]- Note that the CP violation effect reduces the Z — Z — Hi coupling because 
Hi can have an admixture of the CP-odd Higgs state a. The vertical axis of Fig. 4 of |3H] 
describes the Z — Z — h coupling in the case of the MSSM with real parameters, which 
is reduced by a factor sin^(/? — a) in comparison to the SM. The CP violating effects 

^Note that On ~ — sin a, O21 ~ cos a and Hi ^ h for m//± ^ mz- 
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can easily be included by using (On cos/? + 02isin/?)^ instead of sin^(/3 — a). For the 
calculation of the b ^ s'-f width in condition (iii) we use the formula of jlH] including the 
0{as) corrections as given in 120] • (iv) constrains fj, and tan/3 (in the squark sector), (v) 
is the approximate necessary condition for the tree- level vacuum stability ||51j. 

Inspired by the gaugino mass unification we take \Mi\ = 5/3 tan^ ^^^'^^2 and rrig = 
{as{mg)/a2)M2 with rrig = M3. In the numerical study for ti^2 decays we take tan/3,M2, 
m^-j , 771^2 , m^^ , \At\, \Ab\, v?At, V^a,,, V^^' V'u(i) and m^i as input parameters, where mj^ ^ 
and m^^ ^ are the on-shell squark masses. From these input parameters we first calculate 
Mq and Mjj according to the formulae 



Mi 



= ^ + m| ± (m| - m?J2 _ 4^2 _ ^* ^|2^ 
1 2 

— (- — - sin^ 9w) cos 2/5 m| — , (43) 

^0 = 1 ("^l + "^1 T (m| - m|)2 - |^ - fi* cot/5|'^ 
2 

- - sin^ 9w cos 2/5 m\ - ml . (44) 

We resolve the sign ambiguity by assuming either Mq > Mq- or Mq < M^: upper (lower) 
signs correspond to the case Mq > Mfy (Mq < Mq)? With Eq. (jH)) this uniquely fixes the 
mixing angle 9i. Next we calculate M^, using Mq and m^^ and then m^^ and the mixing 
angle 9i as well as the mass eigenvalues and mixing matrices of the charginos, neutralinos 
and the neutral Higgs bosons. For 61.2 decays we take the same input parameters with 
mi^ replaced by m^^ and proceed in the analogous way by interchanging ^ M^. 



3.1 Top squark decays 

In this subsection we present numerical results for the dependence of the ti and ^2 par- 
tial decay widths on ipAti V^Ab and y?^- In order not to vary too many parameters we 
fix (m4^,mt,,m^J = (350,700,170) GeV [(350,800,170) GeV] in the plots for the ti[t2] 
decays. We have selected the parameters in this subsection such that fermionic as well 
as bosonic decays are allowed at the same time. In particular, the choice m^^ = 170 GeV 
has been made to allow the decays ii — >• hiW^ and ti — >• biH^. We consider the cases 
Mq > Mjj and Mq < Mfj, calculating the values of Mq, Mq and corresponding to 
mf^, mi^ and m^^ for each case, as explained above. 

We show in Fig. ^ the partial decay widths and branching ratios for ti — > xt^i ti ~^ 
xt^i ^1 ~^ Xit and ti —>■ W~^bi as a function of (pAt for the parameters tan/5 = 6, 
M2 = 300 GeV, \Ab\ = \At\ = 800 GeV, cp^ = vr, (^u(i) = Va, = 0, m^i = 900 GeV and 
two values of \fi\ = 250 and 350 GeV. Figs. □(a) - (c) (Figs.[I](d) - (f)) are for Mq > Mq 
{Mq < Mq). We first discuss Fig. □ (a), (b) for the case Mq > Mq and = 350 GeV. 
As can be seen in Fig. [T] (a), T(ti Xib) and r(ti x°t) show quite a significant ipAt 
dependence. The corresponding branching ratios are shown in Fig.^(b). For (fAt ~ and 

^The hierarchy is determined by the ii mixing angle 6j:, which can be determined by cross section 
measurements with polarized beams |S1|5]. 
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27r the decay ti Xit dominates, whereas for (fAt ~ the decay ti — > Xi^ has the largest 
branching ratio. This decay pattern can be explained in the following way: For Mq > Mfj 

the ii is tij-like. For > M2 and the parameters chosen, the chargino {xi) is W^^-like 
with m-± = 279 GeV, so that the decay ti — >• xt^ is suppressed by the vanishing tR-h-W^ 
coupling and by small phase space. For the parameters chosen we have \At\ 3> tan/3, 
therefore and hence 61 depend only weakly on (fAf However, we have (pi ~ ipAt (see 

Eq. (0)), therefore T(ti xtb) behaves like (1 — cosifAt)'- the leading coupling term in 
this decay is = — e"""^' cos 6^41^11 + sinO^YtVu (Eq. (fT^ ). which consists of two terms of 
comparable size, the phase (pi{^ v^aJ entering only in one of the two terms. T(ti xtb) 
is very small for cpAt = and 2n because the two terms nearly cancel each other. The ipAt 
dependence of T{ii x^t) is less pronounced compared to T{ii — >• xt^) due to a more 
complex coupling structure (Eqs. (fT^ and (fTTj) ). For this reason B{ti — >• x^t) dominates 
for ipAt ^ OAn and ipAt ^ l-Gn, whereas B{ii xtb) is larger for OAn < (fAt ^ I.Ott. The 
branching ratio of ti W^bi is strongly suppressed for this set of parameters with rather 
small tan (3 = 6 for which bi is almost purely fei^-like. In Fig. ^ (c) we plot the branching 
ratios for \u\ = 250 GeV. In this case the lighter chargino has a mass m-± = 230 GeV and 

a significant higgsino component. Hence the decay ti xtb has a large phase space and 
large amphtude (due to the large top Yukawa coupling gYt) and dominates independently 
of (pAt, resulting in a weak ipAt dependence of the branching ratios. For = 250 GeV 
also the decay channel ti xt^ ^^xt ~ ^"^^ GeV) is open. 

Figs, m (d) and (e) show the partial decay widths and branching ratios of ti — > xtb^ 
xtb) Xit and W'^bi against ipAt foi' < 1/^1 = 250 GeV and the other parameters as 
above. In this case ti is ti-like, therefore for = 250 GeV (see Fig.[T](d)) r(ti — > xtb) is 
about three times as large as for Mq > Mq and = 350 GeV (Fig. [T] (a)). T(ti xtb) 
behaves like (1 — cos (p At), which is again caused by an interplay of the two terms in the 
leading coupling (Eq. (fT^ ). For ipAt ~ the decay ii xtb is suppressed and the 
branching ratios of ti xtb, ii W^bi and ii x^t reach 25%, 22% and 11%, 
respectively (Fig. [T] (e)). For 0.27r < pAt ^ l-Svr the partial decay width and hence the 
branching ratio of ti — > xtb is clearly largest. B(ti — *• x^t) has values around 10%. 
B(ti — » W'^bi) is rather small because bi ~ bn in this case, ti xtb is suppressed by a 
small phase space. In Fig.^(f) we show the corresponding branching ratios for Mq < Mq 
and = 350 GeV. In this case the mixing in the bottom squark sector increases and 
B(ti — >• W^bi) reaches values around 10% even for tpAt ~ ^- The decay ti — >• xtb has 
the largest branching ratio because ti is t^-like and xt i^ almost H^"'"-like. Hence in this 
scenario all branching ratios show a less pronounced phase dependence. In the scenarios 
of Fig. [T] we have calculated also the 9?u(i) dependence of the partial decay widths and 
branching ratios. By inspecting Eqs. (fTBj) - (fTHjl one can see that only r(ti xtt) could 
be sensitive to pv(i)- However, for tan/5 = 6 the p\j(i) dependence is already rather small. 
This results in a weak <y2u(i) dependence of the branching ratios. 

In Fig. 121 we show the tan/5 dependence of B{ti x^t) for M2 = 300 GeV, = 
300 GeV, \Ab\ = |A| = 600 GeV, ¥?u(i) = ^a^ = 0, mH± = 500 GeV, and p^ = 
0, 7r/2, 57r/8, vr with (a) (fAt = and (b) (fAt = assuming Mq > Mq. As can be seen 
this branching ratio is insensitive to (/?^ for tan/5 > 15. This is mainly due to the /i/ tan/5 
dependence of the ti-tR mixing term and the insensitivity of the masses and mixing of 
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r/GeV (a) Mq > M^, = 350 GeV r/ GeV (d) < M^, = 250 GeV 




Figure 1: (a), (d) Partial decay widths F and (b), (c), (e), (f) branching ratios B of 
the decays ti xt^ (sohd), ti — >■ (dotted), ti Xi^ (dashed) and ti — > W^bi 
(dashdotted) for tan/? = 6, Ma = 300 GeV, \Ab\ = |A| = 800 GeV, = tt, (^u(i) = 
cp^^ = 0, rui^ = 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV and mH± = 900 GeV. In (a), 
(b) and (f) the decay ii — > xt^ is kinematically forbidden. 
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Xit) for (p^ = (solid), 7i/2 (dashed), bn/S (dash- 
(a) and vr (b), Ms = 300 GeV, = 300 GeV, 



Figure 2: Branching ratio B{ti 
dotted) and vr (dotted) with ip^t 
\Ab\ = \At\ = 600 GeV, (^u(i) = ^a, = 0, m^-^ = 350 GeV, m^-^ = 700 GeV, m^^ = 170 GeV 
and mj7± = 500 GeV, assuming Mq > Mq. In (a) the case = is excluded by the 
limit B{h — > S7) < 4.5 x 10~^, and the lines for ip^^ = 7t/2 and Lp^ = tt end in full circles 
beyond which B{h 57) > 4.5 x 10"^ for (tan/3 > 21) and B{h 57) < 2.0 x 10"^ for 
(tan/5 > 13), respectively. 



Xi to (p^ for large tan (3. Two curves in Fig. Efa) end in full circles beyond which the 
experimental constraint from B{h — > 57) is violated: in case y?^ = 7r/2 {(p^ = it), one has 
B{b ^ 57) > 4.5 X 10"^ {B{b 57) < 2.0 x 10"^) for tan/5 > 21 (tan/5 > 13). The case 
ip^ = is completely excluded for this set of parameters. However, for (pAt = tt (Fig. El 
(b)) the constraints from B{b 57) are always fulfilled. 

We have also calculated the tan j3 dependence of the branching ratios of the ti decays 
for Mq < Mjj. B{ti — * Xit) is smaller in this case. Therefore the effect of the phase on 
the tan (3 dependence is also smaller than in Fig. |21 Moreover, for Mq < Mjy the situation 
is different from that shown in Fig. |2l because now for ipAt = the whole tan /5 range is 
allowed, whereas for (pAt = tt the constraints from B{b 57) limit the tan/? range. 

In Fig. El (a) we show a contour plot for B(ti Xit) ^ function of ipAt and for 
tan/3 = 6, M2 = 300 GeV, |/i| = 500 GeV, |A| = |A| = 800 GeV, y^uci) = ^a, = and 
mH± = 600 GeV, assuming Mq > Mq. For the parameters chosen the (p>At dependence is 
stronger than the (p>^ dependence. The reason is that these phase dependences are caused 
mainly by the ii -Ir mixing term (Eq. (jH)), where the tp^ dependence is suppressed by 
cot (5. The ip>^ dependence is somewhat more pronounced for ip>At ~ ^ than for ip>At ~ 0, 27r. 
In Fig. El (b) we show the contour plot of B{t\ —>■ Xit) as a function of (pAt and \At\ for 
(p^ = and \At\ = \Ab\. Clearly, the (pAt dependence is stronger for larger values of \At\. 
For Mq < Mjj we have obtained a similar behavior. Note that the phase dependences of 

the decay branching ratios of ti — > xt^, ti xt^ ^'^^ ^1 ~^ Xit analysed in Figs. HI El 
and El (where the decay ti — biH~^ is kinematically forbidden) would be present also for 
bi masses significantly larger than 170 GeV. 

In Fig. Hlwe show the contour plot for B(ti Xi^) as a function of (p>At and ipA^ for 
tan p = 30, M2 = 300 GeV, = 300 GeV, \Ab\ = | A| = 600 GeV, = vr, v?u(i) = and 
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V'm/^ (a) \At\/GeY (b) 




Figure 3: Contours of B{ii Xit) for tan/? = 6, M2 = 300 GeV, = 500 GeV, 
= ^A, = 0, rrii^ = 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV, m^i = 600 GeV, 
with (a) \At\ = \Ab\ = 800 GeV and (b) = 0, |A| = assuming Mq > M^. The 
shaded area marks the region excluded by the Higgs search at LEP (i.e. by the condition 

(ii))- 

= 160 GeV, assuming Mq > Mq. As can be seen, there is a remarkable correlation 
between ipAt and ipAt, which turns out to be relatively independent of y?^. The ipAt-y^A,, 
correlation can be explained by the behavior of the partial decay width T{ti H^bi), 
which influences all decay branching ratios. As bi ~ in this case, the relevant coupling 
for ti H'^bi is Cf- ~ (7?.*G')^2 (see Eq. dHj)). 7?.* depends on (p^t via tL-tn mixing, 

^1^1 

whereas G depends on (fA^ via the coupling term nib {Al tan f3 + n). As ~ <fAt in this 
case, we have (7^*G)l2 ^ e'^'^^^-^^b) cos 6j:-mb\Ah\ tan/3+sin 6^-2mtmb/ sin 2/5 which clearly 
shows the correlation between if At and (fA^ apart from the much weaker ipAt dependence 
of 6^. Note that here the small value for the bi mass (m^^ = 170 GeV) is important: for 
a larger bi mass the decay ti biH^ would not be allowed kinematically and hence the 
ipAf, dependence shown in Fig. 0] would disappear. 

For Mq < Mjj the decay ii — > H^bi dominates for all ipAt and ipA^ resulting in a 
weaker phase dependence of all branching ratios. Hence also the correlation between (pAt 
and (fAt in B(ti xt^) is less pronounced. However, in the scenario of Fig. |3]one has 
B{b S7) > 4.7 X 10-^ for Mq < M^. 

For the heavier top squark ^2 more decay channels are open. Besides the fermionic 
decay modes ^2 X^b, (j = 1, 2; i = 1, . . . , 4) there are also the bosonic decay modes 
^2 — >■ W~^bj, H~^bj, Zii, Hiii (j = 1,2; i = 1,2,3). In Fig. (a) we show the branching 
ratios for ^2 Xi'2^ ^2 — * X2 3 4^ ^ function of ipAt for tan /? = 6, M2 = 300 GeV, 
\^\ = 500 GeV, j^bl = \At\ = 500 GeV, cp^ = (^u(i) = <^a, = 0, m^^ = 350 GeV, 
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Figure 4: Contours of B{ti — > Xi^) as a function of ipAt and ipA^, for tan/5 = 30, M2 = 
300 GeV, = 300 GeV, \Ab\ = \At\ = 600 GeV, iff, = vr, y^uci) = 0, m^-^ = 350 GeV, 
m^^ = 700 GeV, m^^ = 170 GeV and mH± = 160 GeV, assuming Mq > Mfj. The shaded 
areas are excluded by the experimental limit B{h — > 57) > 2.0 x 10~^. 

mi^ = 800 GeV, m^^ = 170 GeV and mH± = 350 GeV, assuming Mq > Mfj. The ipAt 
dependence of B{t2 is again due to a direct phase effect, because the leading 

coupling £2^5 j = I5 2 (Eq. (fT^ ) consists of two terms, with the phase (pi{^ V'aJ entering 
only the factor 71^2 the second term. Therefore, the shape of -6(^2 is like 

(1 ± cos if At)- Also the phase dependence of the branching ratios into neutralinos is 
mainly due to a direct phase effect. In r(t2 — > X?^); = 2, 3, 4 the phase v^f (~ v^aJ enters 
into the second term of the couplings a^^ and b^^ (see Eq. (HHI)) via 7^22- ^(^2 X2^) 
the coupling 032 dominates and the size of its second term is smaller than 10 % of its first 
term. Hence the |a22p term in the width of Eq. 1)211) creates its weak phase dependence like 
(10 + cosy) At)- However, for r(t2 — X^t) the mixing phase enters mainly into the second 
term of the partial width via Re(a23623) ~ ^^(J^22'^2i) ~ cosip^ ~ cos (p At, resulting in 
a shape like (1 + cosv^aJ- For T{t2 x^t) the two terms in have comparable size 
resulting in a strong (fAt dependence of the terms {a^^l"^ and Re(a24624) in the partial 
width which eventually causes the branching ratio to behave like (1 — cos v?At)- 

In Fig. El (b) we show the branching ratios for the bosonic decays ^2 Zti and ^2 
Hiti {i = 1, 2, 3) for the same parameter values as above. The shape of B{t2 Zti) is like 
(1 — cosifAt) 5 which is solely due to the factor | sin26'jp (see Eq. ^^)- Quite generally, 
the phase dependence of r(t2 — > Hkti) is the result of a complicated interplay among 
the phase dependences of the Hk masses, the top squark mixing matrix elements 7^*^-, 
the neutral Higgs mixing matrix elements Oij and the direct top squark-Higgs couplings 
of tLtij0i,2 and tLtjia. In the present example the 9? At dependence of the partial widths 
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Figure 5: if At dependence of branching ratios of the decays (a) t2 (sohd, 
black/gray), X2/3/J- (dashed, black/gray/hght gray) and (b) t2 — > Zti (dashdotdot- 
ted), ^2 — (long dashed, black/gray/hght gray) for tan/3 = 6, M2 = 300 GeV, 
= 500 GeV, \Ab\ = \At\ = 500 GeV, (/?^ = y^uci) = ^a, = 0, m^-^ = 350 GeV, 
rrii^ = 800 GeV, m^^ = 170 GeV and m^i = 350 GeV, assuming Mq > M^. Only the 
decay modes with B > 1% are shown. The shaded areas mark the region excluded by 
the experimental limit B{b 57) < 4.5 x 10^^. 

r(t2 Hi 2;iii) is mainly due to the LfAt dependence of the factors TiJ and Cit^j^Hiiii) in 
Eqs. (I30|) - (jnH), whereas the ipAt dependence of the Oij is less pronounced in this case^. 

We have also calculated the branching ratios of the ^2 decays for Mq < Mjy. In this 
case no constraints on the ipAt range from the B{b ^ 57) data arise in the given scenario. 
The ifiAt dependence of B{t2 — > Zti) and B{t2 — > -f/^1,2,3^1) is very similar to that shown 
in Fig. |S1 (b). The leading branching ratios are now B{t2 xt^)> B{t2 H^h2) and 
B{i2 Vr+62) with the values 17%, 15% and 13% for = 0,27r and 18%, 7% and 
24% for if At = ^5 respectively. 

Furthermore we have calculated the V5u(i) dependence of the branching ratios of the ^2 
for the scenario of Fig. El It turns out to be very weak because r(t2 Xit) (with Xi ~ B) 
is suppressed in this scenario. The (pu(i) dependence stems only from that of r(t2 X?^); 
i = 2,3,4. X2 ^^d X34 are wino- and higgsino-dominated, respectively. Hence the masses 
rriy^o ^ ^ and mixings Nij {i = 2, 3, 4) of X2,3,4 rather insensitive to the bino phase (pu(i). 

3.2 Bottom squark decays 

In the discussion of 61^2 decays we fix tan f3 = 30 because for small tan f3 the bottom squark 
mixing is too small to be phenomenologically interesting. We fix the other parameters as 
m^^ = 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV, mH± = 150 GeV and M2 = 200 GeV. 
We have chosen a relatively small value for the ti mass to allow for the decay bi H^ti 
which has a rather strong dependence on ipA^- 

^For completeness we remark that the effect of the phase dependence of ii — ij — Hk coupHngs also 
shows up in processes like e^e~ — > titiHi |52| . 



15 




Figure 6: y^Ab dependences of (a) partial widths and (b) branching ratios of the decays 
bi x^b (sohd), bi — > x^b (dashed), bi H~ii (dotted) and 61 W^U (dashdotted) 
for tan/3 = 30, M2 = 200 GeV, = 300 GeV, |A| = \At\ = 600 GeV, ip^ = vr, 
= (^u(i) = 0, m^^ = 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV and mH± = 150 GeV, 
assuming Mq > Mj^. 



In Fig. ini we show the partial decay widths and the branching ratios of bi Xi 2^) 
H^ii, W^ii as a function of v^a^ for = 300 GeV, \A}y\ = \At\ = 600 GeV, ip^ = n 
and (pAt = V'u(i) = 0, assuming Mq > M^. In the region 0.57r < ipA^ < l-5vr the 
decay bi — > H^ii dominates. The ipAt dependence of r(6i — > H^ti) is due to the term 
mft|Aft|e~*'^'^!' tan/? in Eq. (|^. The partial decay widths r(6i Xi2^) almost ipA,, 
independent because the (fA^ dependence of the 6-mixing matrix TZ'' nearly vanishes for 
tan P = 30. Hence the ipAt dependence of the branching ratios B{bi Xi 2^) is caused by 
that of the total decay width. r(6i — * W~ii) is suppressed because bi ~ bji and ti ~ 
in this scenario (since also Mq > Mjy). For the scenario of Fig. El the case Mq < Mjj is 
excluded by the experimental lower limit B{b — > 57) > 2.0 x 10~^. 

The V5u(i) dependence of the partial decay widths and branching ratios is very weak 
in the scenario of Fig. [Ul V^u(i) enters only into r(6i x^b) {k = 1, 2) which are nearly 
independent of V2u(i) because 61 ~ bn and hence mainly /i^j^ in a^^ and f^f. in b^j^ contribute 
(see Eqs. (fT^ - (0)- Then the phase of Nki, which strongly depends on 'fiu(i), almost 
drops out in Eq. (PT|) . Furthermore, the masses m^o and mixing matrix elements Nij of 
the Xi'-sector are insensitive to (pu(i) for large tan /?. 

For large tan/3 one expects also a significant \Ah\ dependence of r(6i H^ii) (see 
Eq. ()28p). This can be seen in Fig.[7|(a) where we show the contour plot of B{bi H^ti) 
as a function of \Ah\ and ipA^, for = 300 GeV, = vr, ipAt = fv(i) = and \At\ = \Ab\, 
assuming Mq > Mfj. The ipA^ dependence is stronger for larger values of \Ab\. Although 
Fig. 13 (a) is similar to Fig. |21(b), the 1^4^! and ipA^ dependence in Fig. 13 (a) is now caused 
by the coupling rUb {Al tan /3 + /x) in Eq. (j2HI)- 

In the case Mq < Mjy, we have (&i,ti) ~ (&l,^/?) (since also Mq > Mjy) and hence 
Cfr ~ mt{At cot P + ^*) (see Eq. (|2Zj))- Therefore B{bi H^ii) is nearly independent 
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Figure 7: Contours of B(hi R-ii) for tan/? = 30, Ma = 200 GeV, = 300 GeV, 
y,^ = vr, ¥Pu(i) = 0, m^^ = 350 GeV, m^^ = 700 GeV, m^-^ = 170 GeV, m^i = 150 GeV 
and (a) \Ah\ = \At\, ipAt = and (b) \Ah\ = \At\ = 600 GeV, assuming Mq > Mfj. The 
shaded areas in (a) and (b) mark the regions excluded by the Higgs search at LEP (i.e. 
by the condition (ii)) and by the experimental limit B{b — > s'j) > 2.0 x 10"'^, respectively. 

of (fAi, which leads to contour hues approximately parallel to the ipA^-di-xis. In this case, 
however, nearly the whole parameter space (i.e. the region with ^ 800 GeV) shown 
in Fig. □(a) is excluded by the limit B{b ^ 57) > 2.0 x 10""^. 

In Fig. [7| (b) we show the contours of B{bi — > H^ii) as a function of (fA^ and ipAt foi' 
I At I = \Ab\ = 600 GeV and the other parameters (except (pAt) as in Fig. [7| (a). As can 
be seen, the ipA^-^PAt correlation is even stronger than in Fig. |3] although it has the same 
origin as that of B{ti xt^)- Note that in the given scenario with mH± = 150 GeV the 
constraint on B{b ^ 57) is only fulfilled for a limited range of (pAf The case Mq < Mfj 
is excluded because B[b —>■ sj) is smaller than 2.0 x 10~^ for this case. Moreover we 
want to remark that even for small tan f3 the 61^2 decay branching ratios can be somewhat 
sensitive to fAn and ip^ |38j. 

In case of the 62 decays more decay channels are open. In Fig.|Slwe show the branching 
ratios for the bosonic decays 62 — * W^^^i,2, Zbi, H^ii^2 and -^1,2,3^1 as a function of ipA^ for 

= 350 GeV, \Ab\ = \At\ = 600 GeV, y?^ = ipAt = and (pv{i) = 0, assuming Mq < Mjj. 
The branching ratios of the fermionic decays are nearly independent of (pA^ in this scenario. 
The phase dependence of r(62 — ^ ^^"^1,2) and r(62 — ^ Zbi) is caused solely by the phase 
dependence of the squark mixing angles 6^ and 61 which is very weak in this scenario. The 
strong ipA^ dependence of r(62 H^ii,2) is caused by the term nib (A^tan/5 + /x) in the 
coupling C^^ (Eqs. and (jlHl))- As 62 ~ and ipi ~ (fAt = vr in this case, the dominat- 
ing term in the coupling CS is (Jl^G)i2 — —e~^'^-'^b cos6imb\Ai,\ tan /9 + 2 sin 6'tmtmb/ sin 2/3 
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Figure 8: ifA^, dependences of the branching ratios of the bosonic decays (a) 62 — ^ 
W~ii/2 (dashdotted, black/gray), 62 — ^ H~ii/2 (sohd, black/gray) and (b) 62 Zhi 
(dashdotdotted), 62 -^^1/2/3^1 (long dashed, black/gray/hght gray) for tan/3 = 30, 
M2 = 200 GeV, = 350 GeV, |A| = \At\ = 600 GeV, ip^ = ipA, = vr, (^u(i) = 0, 
m^^ = 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV and mH± = 150 GeV, assuming 



for 62 — ^ H~ii and {Tl^G)2i — —e^^'^^b sin 6'(-m;,|A;,| tan/? — 2cos9fmtmb/ sin 2/? for 62 
H~t2. Therefore, -8(62 H^ti) and -8(62 — ^ H~t2) behave hke (1 + cosyj^ij and 
(1 — cosi/J^J, respectively. As in the example for the ^2 decays (Fig. E)) the v^Aj, depen- 
dence of B{h2 Hihi) (i = 1, 2, 3) is mainly due to the phase factors explicitely appearing 
in Eq. (p?7j) whereas the ipA^, dependence of the Oij is less pronounced. Furthermore, there 
is only a small mixing in the bottom squark sector with 62 ~ &_r and fei ~ 6^ in this 
scenario. Hence the phase dependence of -B(&2 ~^ Hihi) can be explained by the phase 
dependence of Cih^Hjhn) (Eq. (jHTj) ). It turns out that Hi and are nearly CP-even 
Higgs bosons (0i,2) with O^i ~ and \ii\02i ~ |v4;,|0ii {i = 1,3), which results in the 
pronounced ipA,, dependence of -8(62 — ^ -f^i.a&i)- H2 is mainly a CP-odd Higgs boson (a) 
with O12 ~ O22 ~ and sin/?|A5| ^ cos/?|/i|, resulting in the weak ipA^ dependence of 

We have analyzed the 62 decay branching ratios also for Mq > M^,. The ipAt depen- 
dence of 5(62 Zhi) and -B(&2 ~^ -^^1,2,3^1) are similar to those in Fig. |H1 (b), but they 
are smaller by a factor of ~ 3. The other branching ratios are nearly independent of (pA,,- 
However, for the scenario of Fig. |Slthe case of Mq > Mj^ is excluded because B{h 57) 
is smaller than 2.0 x 10~^ for this case. 

The </2u(i) dependence of the partial decay widths and branching ratios in the scenario 
of Fig. IHl with Mq < Mf) is very weak for the same reason as in the scenario of Fig. IHl 

with Mq > Ml) for the decays of the bi. 



4 Parameter Determination 

We now study to which extent one can extract the underlying parameters from measured 
masses, branching ratios and cross sections. Having in mind that the squark masses are 
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relatively large in the scenarios considered, we assume the following situations: (i) A high 
luminosity linear collider like TESLA can measure the masses of charginos, neutralinos 
and the lightest neutral Higgs boson with high accuracy |M1 155j. In the case that the 
squarks and the heavier Higgs bosons have masses below 500 GeV, their masses can be 
measured with an error of 1% and 1.5 GeV, respectively, (ii) For SUSY particles with 
masses larger than 500 GeV their masses can be measured at a 2 TeV e^e~ collider, such 
as CLIC. The masses of heavy Higgs bosons and squarks can be measured with an error of 
1% and 3%, respectively [221 EZl- For the production we can get an e~ beam polarization 
of P_ = 0.8 and an beam polarization of P+ = 0.4. (iii) The gluino mass can be 
measured at the LHC with an error of 3% |53jj. (iv) mt can be measured with an error of 
0.1 GeV. In this case this error can be neglected in the fitting procedure |56 . We assume 
that the error on can also be neglected, (v) The branching ratio of 6 — S7 can be 
measured within an error of 0.4 x 10~^. 

We do not take into account additional information from LHC about the tj and hi 
systems, because the amount of information available strongly depends on the scenario 
realized in nature [SH]. For example in the SPSla scenario the decay channel ti ^ h + xt 
cannot be identified experimentally, because the chargino decays into a scalar tau to 
practically 100%. and hi production at LHC will probably not give enough information 
about the stop and sbottom mixing angles. Moreover, the formulae for the production 
cross sections at LHC, which exist in the literature are for the real case only and do not 
include complex phases, which might be important for the one-loop corrections. This is 
the main reason why we did not consider LHC data for the stop and sbottom systems, but 
data from CLIC. To clarify the situation at LHC concerning the scenarios we considered 
would require additional theoretical work including complex phases and further Monte 
Carlo studies, which are beyond the scope of this paper. 

Our strategy for the parameter determination is as follows: 

(i) Take a specific set of values of the underlying MSSM parameters. 

(ii) Calculate the masses of t,, 6,, Xki ^^e production cross sections for e^e~ 
iiij, and e+e~ hihj, and the branching ratios of the ii and hi decays. 

(iii) Regard these calculated values as real experimental data with definite errors. 

(iv) Determine the underlying MSSM parameters and their errors from the "experimen- 
tal data" by a fit using the program MINUIT ^59j. 

We consider two scenarios in the following, one with small tan/? and one with large 
tan/?. The small tan/? scenario is characterized by: = 169.6 GeV, = 408.8 GeV, 
Mq = 623.0 GeV, \At\ = \Ab\ = 800 GeV, ifA^ = ^a, = vr/4, ipu(i) = 0, M2 = 300 GeV, 
/i = —350 GeV, tan /? = 6, nig = 1000 GeV, and mH± = 900 GeV. (Here we do not assume 
the unification relation between rrig and M2.) The resulting masses and their assumed 
experimental errors are: m-± = (278.5 ± 0.2) GeV, m-± = (384.5 ± 0.3) GeV, m^o = 
(148.7 ± 0.3) GeV, m^o = (277.8 ± 0.5) GeV, m^o = (359.1 ± 0.3) GeV, m^o = (382.0 ± 
0.7) GeV, ruH, = (115.47±0.05) GeV, m^^ = (896!5±9.0) GeV, niH:, = (897.1±9.0) GeV, 
rui^ = (350.0 ±3.5) GeV, m^^ = (700.0 ±21.0) GeV, m^^ = (170.0 ± 1.7) GeV, and m^^ = 
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Table 1: Decay branching ratios (in %) for top squarks and bottom squarks in the two 
considered scenarios. Corresponding values of the underlying MSSM parameters are given 
in the text. 





scenario with tan (3 = 6 


scenario with tan /5 = 30 


channel 




^2 


hi 


&2 


ii 


^2 


h 




qXi 


66.4 


1.6 


100 


0.6 





0.6 


63.5 


0.6 


QX2 





7.5 





8.7 





8.5 


36.1 


10.3 







13.1 





0.3 





11.1 





4.6 


QXa 





6.6 





2.4 





8.7 





4.6 


q'Xi 


33.1 


19.2 





9.7 


100 


22.5 





14.1 


I'xt 





1.6 





21.0 





6.8 





24.2 




0.5 


0.3 





56.8 





3.1 


0.4 


27.1 


H^q'i 

















7.7 





6.4 


Zqi 




26.9 




0.2 




13.1 




1.5 


Hiqi 




23.4 




0.2 




12.7 




1.4 


H2qi 














2.8 




2.7 


H^qi 














2.4 




2.7 



(626.0 ± 19.0) GeV. Moreover, we find B(h — » 57) = 3.6 x 10 ^. The corresponding top 
squark and bottom squark branching ratios are given in Tab. The large tan (3 scenario 
is specified by = 360.0 GeV, Mjj = 198.2 GeV, Mq = 691.9 GeV, \At\ = 600 GeV, 
y,^^ = 7r/4, \Ab\ = 1000 GeV, ^a, = 37r/2, ipu(i) = 0, M2 = 200 GeV, fx = -350 GeV, 
tan (3 = 30, rrig = 1000 GeV, and = 350 GeV. The resulting masses and their 

assumed errors are: m-± = (188.2 ± 0.5) GeV, m-± = (374.2 ± 0.9) GeV, m^o = (98.2 ± 
0.6) GeV, m^o = (188.2 ±0.9) GeV, m^o = (358.5 ± 0.9) GeV, m^o = (371.6 ±2.0) GeV, 
ruH, = (113.63 ± 0.05) GeV, = '(340.7 ± 1.5) GeV, m^a = (341.1 ± 1.5) GeV, 

rrii^ = (210.0 ± 2.1) GeV, m^-^ = (729.0 ± 22.0) GeV, m^^ = (350.0 ± 3.5) GeV, and 
= (700.0±21.0) GeV. Moreover, we have B{b — >■ 57) = 4.4 x 10~^. The corresponding 
top squark and bottom squark branching ratios are given in Tab. H We have chosen a 
relatively small 61 mass in the small tan/3 scenario and a relatively small ti mass in the 
large tan (3 scenario. As a result of this in the two scenarios considered the observables in 
the ti and bi sectors are sufficient to determine all ti and hi parameters. 

We have taken the relative errors of chargino and neutralino masses from OEH], which 
we rescale according to our scenario; in case of tan /3 = 30 we have taken into account an 
additional factor of 3 for the errors (relatively to tan (3 = 6) due to the reduced efficiency 
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in case of multi r final states from decays of charginos and neutralinos as indicated by 
the studies in [7j. 

A detailed Monte Carlo study of the ti production e~^e~ titi and the ti decays 
ii ~^ and ii bxt at TESLA (y/s = 500 GeV and C = 500 fb~^) has been performed 
in ^U] for real MSSM parameters. These results cannot directly be used for our error 
analysis, because we consider additional ii and ^2 decays. To the best of our knowledge 
no Monte Carlo studies exist which include all of the ii and ^2 decays considered in our 
analysis. Therefore, we have taken only statistical errors for the production cross sections 
and branching ratios by calculating the corresponding number of events for the decay 
ti — >• X as 

N = 2C \a{iiti) + aiiih)] B{ti X) (45) 

and analogously for ^2, &i and 62 decays. For definiteness we take an integrated luminosity 
£ = 1 ab~^ at a c.m.s. energy = 2 TeV (i.e. at CLIC). We do not take systematic 
experimental errors for the cross sections and branching ratios into account since we are 
not aware of any study considering the systematic errors. Instead we have doubled the 
statistical errors obtained above. The evaluation of the systematic experimental errors 
would require further Monte Carlo studies for a specific linear collider which, however, 
are beyond the scope of our paper. 

For the determination of the squark parameters we have used the information ob- 
tained from the measurement of the squark masses at threshold and the squark pro- 
duction cross sections at ^/s = 2 TeV for two different (e~,e^) beam polarizations 
{P_,P+) = (0.8,-0.4) and (P-,P+) = (-0.8,0.4). Here we have assumed that a to- 
tal effective luminosity of 1 ab~^ is available for each choice of polarization. The cross 
section measurements are important for the determination of | cos^^p and | cos^^p as can 
be seen from Eq. ()25|1 and the formulas for the cross sections in [H]. In the numerical 
evaluation of the squark production cross sections we have included initial state radiation 
according to [60;,. In addition we have used the information from all branching ratios in 
Table ^ with the corresponding statistical errors. These branching ratios together with 
the masses and cross sections form an over-constraining system of observables for the 
underlying parameters M? , M?, M^, Re(A), Im(A), Re(v4fe), Im(Af,), Re(Mi), Im(Mi), 
M2, Re(/i), Im(/i), tan/5, rrig, and tuh^- The latter two enter the formulas for the neu- 
tral Higgs masses and mixing. We determine these parameters and their errors from the 
"experimental data" on these observables by a least-square fit. The results obtained are 
shown in Table El Note that the sign ambiguity for the imaginary parts of the parameters 
is due to the fact that we consider CP-even observables. This ambiguity can in principle 
be resolved by considering appropriate CP-odd observables (as proposed in [201 ^ E2) 
the analysis. As one can see, all parameters except At, can be determined rather precisely, 
tan f5 can be determined with an error of about 3% in both scenarios. The relative error 
of the squark mass parameters squared is in the range of 1% to 2%. At can be measured 
within an error of 2 - 3% independently of tan jS. The reasons for this are: (i) the mixing 
angle in the top squark sector, which can be measured rather precisely using polarized 

beams, depends strongly on At and (ii) At influences strongly the corrections to the 
mass of the lightest Higgs boson. The situation for Ab is considerably worse: in case of 
small tan/3 one gets only an order of magnitude estimate. The reason is that both the 
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Table 2: Extracted parameters from the "experimental data" of the masses, production 
cross sections and decay branching ratios of tj and 6j. The original parameters for each 
scenario are given in the text. 



scenario 


tan j3 = 6 scenario 


tan /? = 30 scenario 


M? 

D 


f2 88 ± 06) X 10^ 


fl 30 ± 02) X 10^ 


M? 

u 


(1.67 ±0.04) X 10^ 


(3.93 ±0.12) X 10^ 


Ml 
Q 


(3.88 ±0.04) X 10^ 


(4.79 ±0.04) X 10^ 


Re(AA 


565.0 ± 13.0 


424.0 ± 14.0 


lm{At) 


±566.0 ± 14.0 


±425.0 ± 15.0 


Re(Ab) 


620.0 ± 190.0 


6.5 ±420.0 


ImfA) 


±230.0 ± 580.0 


±999.0 ± 52.0 


Re(Mi) 


149.3 ± 0.3 


99.6 ± 0.6 


Im(Mi) 


1.0 ± 1.5 


-0.5 ± 2.8 


M2 


300.0 ±0.4 


200.0 ±0.5 


Re(/i) 


-350.0 ±0.3 


-350.0 ±0.6 


Im(yu) 


-0.02 ±0.9 


1.5 ±5.0 


tan /? 


6.0 ±0.2 


30.0 ±0.8 


m-g 


1000.0 ±30 


1000.0 ± 30 


mH± 


900.0 ±5.0 


350.0 ±0.8 



bottom squark mixing angle and the bottom squark couplings depend only weakly on Ah 
for small tan /3. In case of large tan (3 the situation improves somewhat in particular for 
the imaginary part of Ab. The main sources of information on Ah are the branching ratios 
of the decays of the heavier bottom squark into a Higgs boson plus the lighter bottom 
squark because the corresponding couplings depend significantly on Ah (see Eqs. flS^ - 
(I37|l ). From this we conclude that the situation for Ah improves in scenarios where these 
branching ratios are large. An additional source of information could be the polariza- 
tion information of the fermions in bottom squark decays as proposed in ^2]- We have 
found that the analogous fit procedure for real MSSM parameters gives a larger value 
for x^'- = 286.6 for the scenario with tan/3 = 6 and Ax^ = 22.5 for the scenario 

with tan/? = 30. In Table 2 most of the central values of the fitted parameters are the 
same as their input values because we have taken the observables calculated from the 
input parameters as "experimental data" . We have checked that a shift within la of the 
"experimental data" leads to almost no change of the errors of the parameters. 



22 



The results presented in Table El depend clearly on the assumed experimental errors 
which have been summarized in the beginning of this section. It is clear that further 
detailed Monte Carlo studies including experimental cuts and detector simulation are 
necessary to determine more accurately the expected experimental errors of the observ- 
ables for our scenarios, in particular the errors of the top squark and bottom squark decay 
branching ratios. Such a study is, however, beyond the scope of this paper. Furthermore, 
an additional source of uncertainty is the theoretical error due to higher order corrections 
etc. |431 144j . We have not taken into account these effects because most formulae given 
in the literature are only for real parameters. Instead we have studied how our results for 
the errors of the fundamental parameters are changed when the experimental errors of 
the various observables are changed: we have redone the procedure doubling the errors of 
the masses and/or branching ratios and/or cross sections. We find that the errors of all 
parameters are approximately doubled if all experimental errors are doubled. Moreover, 
in this way we can see to which observables an individual parameter is most sensitive. 
We find that precision on the top squark parameters At, and M?. is sensitive to the 
accuracy of the top squark mass measurement at the threshold as well as to the precision 
of the measurement of the total top squark pair production cross sections in the contin- 
uum using polarized beams. The error of At is also very sensitive to the error of the 
lightest Higgs boson mass due to the large top squark loop corrections. The precision 
on the parameters M'^ and is sensitive to the accuracy of the bottom squark mass 
measurement. The accuracy of Ab is most sensitive to the precision of the measurements 
of the branching ratios for the bottom squark (and top squark) decays into Higgs bosons. 
The precision of /i is more sensitive to the errors of chargino and neutralino masses than 
to the errors of the top squark and bottom squark observables. In the case of large tan 
the precision of tan f3 depends to some extent on the precision of the bottom squark pair 
production cross sections and to a lesser extent also on that of the bottom squark decay 
branching ratios. 

For the determination of the ti and 6j parameters the measurements of the branching 
ratios of the squark decays into Higgs bosons together with those of the squark mixing 
angles from the production cross sections are important. Therefore, we need to obtain 
information about ti, ^2, &i and 62 production and decays separately. This can be achieved 
at a Linear Collider by suitable choices of the c.m.s. energy. We remark that, in the case 
rrit^ , m^^ > 500 GeV the measurements of the cross sections, masses and branching ratios 
of t2 and 62 at an e~^e~ linear collider with ^/s = 2 TeV are necessary for the determination 
of At and Afy, otherwise this might not be possible. However, additional information from 
the LHC on the tj and 6j masses and some of the decay channels would certainly improve 
the situation. In the error estimate presented here we have assumed that many decay 
channels of the ti and bi are open. If this is not the case, then the missing information 
could be obtained by studying the decay properties of the heavier charginos, neutralinos 
and Higgs bosons into ii and bi. 
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5 Summary 



In this paper we have studied the decays of top squarks and bottom squarks 6j in 
the MSSM with complex parameters At, Af,, n and Mi. We have taken into account 
the exphcit CP violation in the Higgs sector induced by ii and bi loops in the case At^i, 
and II are complex. We have presented numerical results for the fermionic and bosonic 
decay branching ratios of ti and bi {i = 1,2). We have analyzed their MSSM parameter 
dependence, in particular the dependence on the CP phases ip^^, cp^ and V?u(i)- We 
have found that the experimental data of the branching ratio of the decay 6 — > 57 can lead 
to considerable restrictions on the MSSM parameter space. In the case of ii decays the 
strong dependence on (f^^ and </7^ is due to the phase dependence of the mixing angle 9^, of 
the mixing phase factor e*"^* and of the Higgs couplings Gu {— - ), G21 (= ~ ) and 

Cit^j^Hitii). In the case of bi decays there can be strong (fAi, dependence if tan/? is large 
and the decays into Higgs bosons are allowed. If the parameters At, Ah, and Mi are 
complex and there is mixing between the CP-even and CP-odd Higgs bosons, the decay 
pattern of ii and bi is even more complicated than that in the case of real parameters. 
This could have important implications for ii and bi searches at future colliders and the 
determination of the underlying MSSM parameters. 

We have also estimated what accuracy can be expected in the determination of the 
underlying MSSM parameters by a global fit of the observables (masses, branching ratios 
and production cross sections) measured at typical linear colliders with polarized beams. 
We have considered two scenarios with tan/3 = 6 and tan/3 = 30. Under favorable 
conditions the fundamental MSSM parameters except At^b can be determined with errors 
of 1 % to 2 %, assuming an integrated luminosity of 1 ab~^. The parameter At can be 
determined within an error of 2 - 3 % whereas the error of Ai, is likely to be of the order 
of 50%. 
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Appendix 



A Masses and mixing in the neutral Higgs sector 



In the complex MSSM the exphcit CP violation in the Higgs sector is mainly induced 
by t and b loops resulting in a 3 x 3 neutral Higgs mass matrix with a mixing of the 
CP-even Higgs bosons 0i and 02 and the CP-odd Higgs boson a. At one-loop level the 
amount of mixing of CP-even and CP-odd Higgs states is approximately proportional to 
sm{ipA_t^ + (Pfj.). The three neutral mass eigenstates are denoted as Hi {i = 1,2,3) with 
masses rriHi < "^^^2 < "^Hs (following the notation of [2H1)- The real orthogonal mixing 
matrix in the neutral Higgs sector is denoted by a 3 x 3 matrix O: 



O 



T 



V ^ / 



(46) 



where 0i, 02 and a are related to the neutral entries of the two Higgs doublet fields by 
H'l = l/\/2{vi + 01 — ^^i), H2 = l/\/2(f2 + 02 + ^0,2) and a = — sin/5ai + cos/?a2. We 
take the parameter ^ = as in j2H]- We have included the full one-loop corrections to the 
Higgs mass eigenvalues m//. and the mixing matrix Oij as implemented in the program 
FeynHiggs2.0.2 [20] • We use these results for mj^. and Oij in our tree-level formula for 
the ^2, &2 decay widths (Eq. (j^ ) and in the constraint (ii). 



B Chargino Masses and Mixing 

At tree-level the chargino mass matrix in the weak basis is given by |^ 

/ M2 V2mwS(^ ^ 



M 



c 



(47) 



C/3 and Sf) are cos j3 and sin /3, respectively. This complex 2x2 matrix is diagonalized by 
the unitary 2x2 matrices U and V: 



U*McV^ = diag(mj±, 



< m-± < m-±. 

Al A2 



(45 



We have neglected one-loop corrections to the chargino mass matrix A4c, as have been 
given in ^44^ IH^ for real parameters. 
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C Neutralino Masses and Mixing 

At tree-level the neutralino mass matrix in the weak basis {B, W^, if °, ii2) is given as 

Binii: 



M 



N 



-mzswcp mzSwSfs 



mzcwcp -mzCwS/3 



\ mzswsp —mzcwsp — |/i|e*'^'' 



(49) 



/ 



where V2u(i) is the phase of Mi, and cw and sw are cosOy/ and smOw, respectively. This 
symmetric complex mass matrix is diagonalized by the unitary 4x4 matrix A^: 



N*MnN'^ = diag(m^o, 



< m^o < 

A.1 



(50) 



We have not included one-loop corrections to the neutralino mass matrix AiN, like those 
given in 15^ for real parameters. 
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